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Abstract 

We investigate the possibility to extract the symmetry energy from multifragmentation data. 
The applicability of the grandcanonical formula earlier proposed by Ono et al. [Phys. Rev. C 68, 
051601(R)] in the case of finite excited nuclei is tested within a microcanonical framework. Rela- 
tively good results are obtained except for large residual nuclei, especially when large sources are 
highly excited. Effects of secondary particle emission and the extent in which relevant information 
may be inferred from experimental observables are finally discussed. 

PACS numbers: 25.70.Pq Multifragment emission and correlations, 24. 10. Pa Thermal and statistical models 
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I. INTRODUCTION 



Isospin dependent phenomena are attracting increasing interest as they hold the promise 
of reveahng the asymmetry term of the nuclear equation of state. At normal nuclear den- 
sity this quantity dictates the structure of the neutron-rich and neutron deficient isotopes, 
while in other domains of density and temperature its behavior is reflected in a variety of 
astrophysical phenomena, such as the structure and evolution of neutron stars and the dy- 
namics of supernovae explosions. In nuclear multifragmentation reactions, the asymmetry 
term influences the neutron-proton composition of the break-up fragments. 

Interpreting multifragmentation in the light of first-order phase transitions in multicom- 
ponent systems, the neutron enrichment of the gas phase with respect to the liquid phase 
comes out as a natural consequence of Gibbs equilibrium criteria and a connection between 
phases chemical composition and the symmetry term can be established Interest- 
ing enough, the phenomenon of isospin fractionation which is systematically observed in 
analyses of multifragmentation data 0, 0, 0, 0, Q], seems to be a generic feature of phase 



separation independent of the equilibrium Gibbs construction j8[ . Indeed, dynamical models 
of heavy ion collisions ^ [lll.ll2[ lisj ] where fragment formation is essentially ruled by the 
out of equilibrium process of spinodal decomposition also exhibit fractionation. Adop ting 
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an equilibrium scenario for the break-up stage of a multifragmenting system, Ono et al 
derive an approximate grandcanonical expression which connects the symmetry term with 
the isotopic composition of fragments obtained in the break-up stage of two sources with 
similar sizes in identical thermodynamical states and differing in their isospin content. 



C. 



sym 
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under the hypothesis that the isotopic distributions are essentially Gaussian and that the 
free energies contain only bulk terms. Here, ai2 is the isoscaling slope parameter, Zi/Ai 
stands for the isospin asymmetry of a fragment produced by the source i{= 1, 2) and T is 
the temperature of the decaying systems. 

In the limit of vanishing temperature, fractionation can be neglected and Zi/Ai can be 
replaced by the corresponding quantity of the sources Q giving 
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of heavy ion col. 
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This expression, which was first proposed in the framework of the EES model Q], as well 
as Eq. has been used extensively on experimental data and results ranging from values 
compatible with the ground state bulk symmetry term 25 MeV) to about half this value 
have been interpreted as sign of dilute matter at freeze-out. Thus, Csym=^4: MeV has been 
obtained in FRS experiments using mid-peripheral collisions of 1 GeV/nucleon ^^"^'^^^Xe 
beams on Pb targets values decreasing from 25 to 15 MeV have resulted from a variety 
isions induced by few tens MeV/nucleon projectiles studied by the Texas 
1^,13,13) while from the fragmentation of excited target residues produced 
in the peripheral collisions of ^^C at 300 and 600 MeV/nucleon incident energies 

the INDRA-ALADIN collaboration extracts a symmetry coefficient which decreases from 25 
MeV for peripheral collisions to 15 MeV for central collisions js^. 

If these values really correspond to the break-up stage of the decay, the implications are 
dramatic. First, if the primary fragments at the equilibrated freeze-out are diluted, sta- 
tistical models which described successfully a wealth of experimental data over more than 
two decades should be completely reformulated. Starting from the geometrical extension of 
fragments which dictates the free volume, to the binding and internal excitation energies 
which enter energy balance and affect fragments partitions, all quantities need serious recon- 
sideration. If, on the other hand, this diluteness would reflect inter-fragment interactions 

n 

at break-up, as advanced by Botvina et al. 7], the situation would be even more severe 
as it would refute the fundamental hypothesis of statistical models, namely the lack of any 
fragment interaction except the Coulomb one. 

Given the implications of such measurements, the present work aims to investigate the 
possibility to extract the symmetry term from multifragmentation data using Eq. (^. Stick- 
ing to the equilibrium hypothesis, Ono's formula has three drawbacks: a) it is grandcanonical 
while the grandcanonical approximation is known to be acceptable only for relatively light 
fragments emitted by large systems with high excitation energies, b) it does not account for 
full mass dependence of the binding energy and c) it holds for the break-up stage of the 
reaction, impossible to access experimentally. Eq. (j21), which was employed to obtain Cgym 
values as low as ~ 10 MeV [20], contains the additional approximation of neglecting isospin 
fractionation. 

A quantitative estimation of the possible deformations induced by these effects can be 
done in the framework of a microcanonical equilibrium model The less serious objection 
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against Eq. (|T]), possible deviations produced by omitting the contributions of surface 
Coulomb and asymmetry terms, may be overcome taking into account the full dependence 
of the binding energy on the system size and will be addressed first. 

Thus, adopting the grandcanonical expression of the isotopic yield of an emitted cluster 
with neutrons and Z protons. 



Y{N, Z) (X exp (S(iV, Z) + /i„iV + /i^Z) 



(3) 



where /i^ and /ip stand for neutron and proton chemical potentials and B{N, Z) represents 
the binding energy, we use as in Ref. the Gaussian shape of Y{N, Z) distributions as a 
function of {Z) when Z (N) is kept fixed to approximate the value of (Z) corresponding 
to the maximum of the distribution with its average value. Under these assumptions, the 
equation 

d{\nY{N, Z)) 
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T 



dB{N,Z)^ 

Z + fJ-n 



(4) 



dN 

applied for two isospin different similar nuclei in identical equilibrium states, relates the 
isoscaling slope parameter a with the difference of partial derivatives of B{A, Z) as a function 
of N. 

If we consider for the binding energy a liquid drop parameterization including surface 
and Coulomb terms, 

B{A, Z) = {ayA - asA^/^) - (a^a^A - 0^0,^^/^)/^ + acZ'^/A^^^ + aaZ^/A, (5) 

and account for full mass dependence of the bulk+surface, isospin-dependent and Coulomb 
contributions, the resulting equation will acquire additional terms corresponding to surface, 
Coulomb and asymmetry energies. 



The < A >~" {n > 0) dependence of these quantities rend the corrections neghgeable 
for heavy fragmentation products and quantitative results will be presented for a case 
relevant to most multifragmentation reactions. 

The approximations of Eq. (Q) and Eq. (0) may be jud ged within a microcanonical 
multifragmentation model 21 1. We use the MMM version j23[ where the equilibrated mul- 
tifragmenting source is characterized by its mass As, charge Zg, excitation energy E, total 
momentum P, total angular momentum L and freeze-out volume V. All configurations 
C : {{Ai, Zi,ei,Pi,Ti}, {i = 1, NF}} permitted by the specific microcanonical conserva- 
tion laws are spanned by a Metropolis trajectory and average values of physical quantities 
are calculated numerically on the basis of the statistical weight of each configuration C. 
Break-up fragments are treated as normal nuclear density malleable objects described by a 
zero-temperature binding energy as in Eq. (jSj) where at,=15.4941 MeV, 0^=17.9439 MeV, 
aj=1.7826, ac=-0.7053 MeV and 0^=1. 1530 MeV |2j]. This liquid drop description is con- 
sistent with a semiclassical Thomas- Fermi approximation [2^! or hot Hartree Fock 
where the effect of temperature is a modified occupation of the single particle eigenstates 
of the mean field hamiltonian. The finite temperature fragment energy functional in this 
approach is thus modified respect to the ground state only for the internal excitation energy 
(e) coming from the occupation of continuum states. To avoid double counting of the free 
particles states the internal energy is cut at the binding level. In our microcanoni- 
cal formalism, the temperature is univocally defined through the thermodynamic equality 
T^^ = dS/dE =< {3N/2 — A)/K >, where K represents the thermal kinetic energy, is the 
product multiplicity and the last equahty stems from the equipartition theorem applied to 



a non interacting system of clusters 



28l |. It is worthwhile mentioning that, due to the high 



energy cut-off in the cluster level densities, the internal fragment excitation energy cannot 
be used to estimate the temperature, and a thermometer based upon internal fragment prop- 
erties would severely underestimate the thermodynamic temperature. In the experimental 
evaluation of Cgym conversely, the temperature is also estimated from fragment properties 
or model calculations, which may lead to an extra source of uncertainty. 

The break-up stage of the decay which contains all information relevant for the equation 
of state is completed with a particle evap oration stage which simulates the subsequent 
disintegration of the excited fragments 23[. This step is important as well as it shows the 
capacity of experimentally measurable quantities to access the physics at break-up. 
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The predictive power of Eqs. (^J and (jH)) has been systematically checked by confronting 
their results with the input symmetry energy which enters fragment definition via the binding 
energy. We have considered a variety of situations in which the size of the emitting sources 
was varied between A = 230 and A = 100, the freeze-out volume covered the usually accepted 
interval {V = 4Vq — 8Vo) and the excitation energy ranged between 2 and 10 MeV/nucleon. 
The lower limit of the source size for which Eqs. (Q), © and (j21) hold was conditioned by the 
possibility to calculate a out of the isotopic composition of light emitted fragments. Indeed, 
for small sources a manifests a relatively strong dependence on the emitted cluster size, its 
calculation as an average value getting disputable [3]. For completeness, the stability of 
the above mentioned formulas was checked against sources isospin variation as well. The 
conclusions are the same and only few illustrative cases will be considered in the following. 

Fig. [T] illustrates the magnitude of the additional terms of Eq. (0) corresponding to 
corrections due to surface. Coulomb and isospin energies as a function of the charge of 
the considered fragment in a case typical for nuclear multifragmentation reactions, (210, 
82) and (190,82) with V = AVq (Vq is the volume at normal nuclear density) and E = 6 
MeV/nucleon. The values of the order of unity rend these corrections negligeable even for 
light fragments, in agreement with expectations of Ref. For this reason, the discussion 
on the possibility to extract the symmetry energy from multifragmentation data will address 
only the original Eq. (^. 

Fig. |2lpresents break-up stage predictions of MMM for the asymmetry term as a function 
of the emitted fragment charge. Results of Eqs. (0) and are shown for the multifrag- 
menting nuclei ((130, 50) and (110, 50)) (upper panel) and ((210, 82) and (190, 82)) (lower 
panel) with V = 4Vq and excitation energy ranging from 2 to 10 MeV/nucleon. The MMM 
results (symbols) are compared to the input symmetry energy of the model, OjO^ — aja^A"^/^ 
(solid curve). As one may notice, Eq. shows a remarkable overall stability against ex- 
citation energy variations but its behavior while modifying the source size or excitation 
energy deserves a more attentive investigation. Firstly, Eq. (0) shows a systematic over- 
estimation of the real value by up to 4 MeV for the largest fragments. With the increase 
of the source size and excitation energy, the overestimation slightly diminishes so that, for 
the Z=82 sources and the highest considered excitation energy, 10 MeV/nucleon, the calcu- 
lated Csym exceeds by less than 1 MeV the real value. This is an interesting manifestation 
of the applicability conditions of grandcanonical approaches in the case of small systems. 
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FIG. 1: (Color online) MMM predictions for the additional surface (dashed), Coulomb (dotted) 
and isospin (dot and dash) contributions to the symmetry energy as a function of Z at break-up. 
The solid curve depicts the sum of the above mentioned terms. The equilibrated systems are (210, 
82) and (190, 82) with V = 4Vo and excitation energy 6 MeV/nucleon. 

Thus, for relatively light fragments emitted by highly excited large systems, grandcanonical 
formulas give reasonable values while more modest results are obtained at low excitation 
energies and for fragments commensurable with the source size. The bending of the curves 
calculated with Eq. (0) for the Z=50 sources at 2 and 4 MeV/nucleon is due to the mass 
and charge conservation specific of a microcanonical ensemble. Similar effect manifests also 
for the Z=80 case, but the limited charge domain in the right panel of Fig. |21 hides it. The 
residual slight overestimation of the input value of Cgym in a domain where the grancanon- 
ical approximation should in principle be acceptable remains an open question and will be 
addressed in the future. However, we have checked that it does not depend on the value 
of the high energy cutoff exp(— e/r) of the fragment internal density, as similar results are 
obtained when r — > oo. 

As expected, Eq. (0) which approximates the isospin content of fragments to the one of 
the source gives good results only for very heavy fragments and low excitation energy. How- 
ever in the fragmentation regime it leads to a systematic underestimation of the symmetry 
energy. In particular. Fig, p s uggests that the low value Csym ~ 15 MeV extracted in some 



experimental analyses 



15l |20| may be compatible with standard ground state values for the 



symmetry energy. In this respect it is interesting to remark that the use of Eq. (jSj) produces 
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FIG. 2: (Color online) MMM predictions for the symmetry energy as a function of Z at break-up. 
The equilibrated systems are ((130, 50) and (110, 50)) (upper panel) and ((210, 82) and (190, 82)) 
(lower panel) with V = 4Vo and excitation energies ranging from 2 to 10 MeV/nucleon. Solid 
curve: input symmetry energy of the model. Cgym calculated according to Eq. (Q) is represented 
with symbols while results of Eq. ((2} are illustrated with dotted and dashed lines. 

an unphvsical apparent reduction of Csym with increasing excitation, similar to the findings 
of Ref. j20|. The explanation of this evolution in the model lays in the steep diminish of a 
not compensated by the temperature increase Q. 

Fig. 121 was obtained for a specific choice of the fragmenting sources. However it was 
proved in Ref. j^l that the isoscaling parameter a manifests a complex dependence on the 
considered pair of sources and their equilibrium state under the microcanonical constraint. 
This fact requires a detailed investigation of the persistence of Eq. under various freeze- 
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FIG. 3: (Color online) MMM predictions for the symmetry energy as a function of Z at break- 
up as obtained using Eq. (open symbols). The stability against freeze-out volume variation 
{y = 4Vb and V = 8Vb) is illustrated for the equilibrated systems (190, 82) and (210, 82) with 6 
MeV/nucleon excitation energy (upper panel). The stability against sources isospin modification is 
checked for three pairs of Z=82 nuclei with V = 4Vo and 6 MeV/nucleon excitation energy (lower 
panel). Results of Eq. @ are indicated using dotted and dashed lines. 

out volumes, source isospin combinations and source sizes. Therefore, Fig. |21 presents the 
results of Eq. ((H) for the same sources ((210, 82) and (190, 82)) with 6 MeV/nucleon 
excitation energy when the freeze-out volume has different values (V^ = 4Vo and V = SVq) 
(upper panels), while the behavior with respect to sources isospin modification is represented 
considering three pairs of Z=82 nuclei ((190,82), (200, 82) and (210, 82)) at V = AVq and 
6 MeV/nucleon excitation energy (lower panels). As in the previous case, the solid curve 
illustrates the ( A ^/^) term. The conclusions confirm the stability of Eq. ^ 

while modifying source isospin and freeze-out volume. A more increased predictive power is 
expected for larger sources with a more advanced fragmentation, where the grandcanonical 
approximation is more reasonable. 
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If the preceding figures suggest that for large systems and high excitation energies, Eqs. 
(0) and © niay allow to extract the symmetry energy coefficient at the fragmentation stage 
from isoscaling observables, it is important to stress that not only the average fragment 
isotopic composition, but also the isoscaling parameter a should be known at the time of 
fragment formation, while the measured a value may have been distorted by sequential decay. 
The effect of secondary decay on the isoscaling parameter is explored in Fig. |3] Since in the 
MMM model the value of a is larger in the asymptotic stage of the reaction than at break- 



29[, Eqs. (0) and © lead to an overestimation of the symmetry energy coefficient. 



up 

If fractionation is neglected using Eq. these two errors tend to compensate and the 
extracted values of Cgym range from 20 to 25 MeV, in better agreement with the real value 
(if sufficiently heavy isotopes are analyzed for the surface influence to be negligible, which is 
not the case in present fragmentation data). It is important to stress that this compensation 
entirely relies on the evolution of a from hot to cold fragments. Unfortunately, in this respect, 

n n 

information in the literature is controversial. Thus, MSU-SMM |^ and IQMD |32|| p lead 
in favor of negligible contribution of sequential evaporation on ot as stated in Refs. [l^ and 
0, respectively. Stochastic mean field Q, AMD 3 and the Markov-chain SMM Q 
wit h C^,jrn. > 15 MeV indicate a values decreasing from the break-up to the asymptotic stage 
QB- Finally, MMM, EES and the Markov-chain SMM with Csym <15 MeV 

Hiict the opposite effect: a values increasing from the break-up to the asymptotic stage 
0,0. The origin of these discrepancies could by due to isospin population of break-up 
fragments, their excitation energy and secondary decay procedure as well, but a pertinent 
discussion of this issue goes much beyond the limits of present paper. 

To summarize, we investigated the possibility of inferring the symmetry energy coefficient 
from isoscaling observables in multifragmentation reactions. Our results indicate that, in 
the limit of large systems and high excitation energies where the use of a grandcanonical 
approach is reasonable, Eq. ((1} is a reliable tool to determine Cgym out of the isotopic 
composition of the break-up fragments but overestimations of few MeV are expected for 
the systems usually produced in multifragmentation reactions. In particular low values of 
Csym respect to the standard saturation density bulk value may be interpreted as an effect 
of surface contributions to the symmetry energy. 

Conversely, if we assume that the density of break-up fragments is close to their ground 
state density [2l|, the symmetry energy coefficient can be considered as known and isoscaling 
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FIG. 4: (Color online) MMM predictions for the symmetry energy as a function of Z in the 
asymptotic stage of the decay as obtained using Eq. (open symbols). The equilibrated systems 
are (190, 82) and (210, 82) with V = 4Vo and excitation energies ranging from 2 to 10 MeV/nucleon. 
Results of Eq. @ are indicated using dotted and dashed lines. In all cases, asymptotic values of 
a and break-up values oi Z/ < A> have been used. 

measurements can be used to probe fractionation. In particular spinodal decomposition is 
predicted to lead to higher fractionation than phase separation at equilibrium jsSi, and this 
effect could be accessed through isoscaling observables. To realize this ambitious program it 
is essential to experimentally control the effect of side feeding on the a isoscaling parameter. 
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